Spontaneous ordering against an external field in nonequilibrium systems 



J. C. Gonzalez- Avella, 1 M. G. Cosenza, 2 V. M. Egui'luz, 1 and M. San Miguel 1 

1 IFISC Instituto de Fisica Interdisciplinar y Sistemas Complejos (CSIC-UIB), E-07122 Palma de Mallorca, Spain 
z Centro de Fisica Fundamental, Universidad de Los Andes, Merida, Merida 5251, Venezuela 

(Dated: November 18, 2008) 

We study the collective behavior of nonequilibrium systems subject to an external field with a 
dynamics characterized by the existence of non-interacting states. Aiming at exploring the generality 
of the results, we consider two types of models according to the nature of their state variables: (i) 
a vector model, where interactions are proportional to the overlap between the states, and (ii) a 
scalar model, where interaction depends on the distance between states. In both cases the system 
displays three phases: two ordered phases, one parallel to the field, and another orthogonal to the 
field; and a disordered phase. The phase space is numerically characterized for each model in a 
fully connected network. By placing the particles on a small-world network, we show that, while a 
regular lattice favors the alignment with the field, the presence of long-range interactions promotes 
the formation of the ordered phase orthogonal to the field. 

PACS numbers: 89.75.Fb, 87.23. Ge, 05.50.+q 



A rather general question considered in the frame- 
work of statistical physics of interacting particles (par- 
ticles, spins, agents) is the competition between local 
particle-particle interactions (collective self-organization) 
and particle interaction with a global externally applied 
field or with a global mean field [H, 0]. Common wis- 
dom answer to this question is that a strong external 
field dominates over local particle-particle interactions 
and orders the system by aligning particles with the bro- 
ken symmetry imposed by the field. However, this is 
essentially an equilibrium concept which is not generally 
valid for generic non-potential interactions. 

In the context of studies of collective phenomena in 
general systems of interacting particles, including biolog- 
ical and social systems, new forms of particle-particle and 
particle-field interactions are being considered. There are 
forms of interaction for which it turns out that a suffi- 
ciently intense external field induces disorder in the sys- 
tem 0, H[, m contrast with the behavior in, for ex- 
ample, Ising-type systems. Other intriguing dynamical 
phenomenon is the collectively ordering in a state dif- 
ferent from the one preferred by the forcing field. The 
external field might break symmetry in a given direction, 
but the system orders, breaking symmetry in a different 
direction. In this paper we examine this situation show- 
ing that these phenomena happen in two recently well 
studied non-equilibrium models [H, H 0, B • What is 
common to these two models is that the particle-particle 
interaction rule is such that no interaction exists for some 
relative values of the states characterizing the particles 
that compose the system. A subsidiary question is the 
dependence of this phenomenon on the topology of the 
network of interactions. We show that the phenomenon is 
not found for particles interacting with its nearest neigh- 
bors in a regular lattice, but it occurs in a globally cou- 
pled system: it emerges as long range links in the network 
arc introduced when going from the regular lattice to a 
random network via small world network (lfj| . 

The vector model, based in the dynamics of cultural 



dissemination of Axelrod model, consists of a set of N 
particles located at the nodes of an interaction network. 
The state of particle i is given by a F-component vector 
C{ (f = 1, 2, . . . , F) where each component can take any 
of q different values C( £ {0, 1, . . . , q — 1}, leading to 
q F equivalent states @. The external field, defined as a 
F-component vector e {0, 1, . . . , q — 1}, can interact 
with any of the particles in the system. 

Starting from a random initial condition, at any given 
time, a randomly selected a particle can either interact 
with the external field or with one of its neighbors. The 
dynamics of the system is defined by iterating the follow- 
ing steps: 

1. Select at random an particle i. 

2. Select the source of interaction: with probability 
B the particle i interacts with the field, while with 
probability (1 — B) it interacts with one of its near- 
est neighbors j. 

3. The overlap between the selected particle and the 
source of interaction is the number of shared com- 
ponents between their respective vector states, d = 

where X f = M f if the source of in- 



C(,Xf> 

teraction is the field, or X? = Cj if i interacts with 
j. If < d < F, with probability d/F, choose h 
randomly such that Gf ^ X h and set Gf = X h . 

The strength of the field is represented by the pa- 
rameter B G [0, 1] that measures the probability for the 
particle-field interactions. In the absence of an external 
field, B = 0, the system reaches a stationary configura- 
tion in any finite network, where for any pair of neighbors 
i and j, d(i,j) = or d(i,j) = F. A domain is a set of 
connected particles with the same state. An homoge- 
neous or ordered phase correspond to d(i,j) = F, 
and obviously there are q F equivalent configurations of 
this state. An inhomogeneous or disordered phase con- 
sist of the coexistence of several domains. In order to 
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FIG. 1: (a) 5 as a function of q for the Axelrod model on a 
fully connected network for B — (solid circles); B — 0.005 
(diamonds); B — 0.05 (empty circles); B = 0.5 (squares); 
5 = 1 (stars). The continuous line is the analytical curve 
1 — (1 — l/q) F , while the dashed line corresponds to the curve 
(1 — l/q) F . (b) a versus q for B = 0.8 (circles) and B = 0.1 
(squares). The values of q c and q* are indicated for B = 0.8. 
Parameter values are N = 2500, F = 10. 



characterize the ordering properties of this system, we 
consider as an order parameter the normalized average 
size of the largest domain S formed in the system. In 
finite networks the dynamics displays a critical point q c 
that separates two phases: an ordered phase (S ~ 1) 
for q < Q r , and a disordered phase (S < 1) for q > q c 

Ida miii. 

First, we analyze the model in a fully connected net- 
work. In the absence of field, i.e. B = 0, the system 
spontaneously reaches a homogeneous state for values 
q < q c w 10 4 (Fig. 1-a ). For B -> and g < q c , the 
field is able to impose this homogeneous state to the 
system, as in a two-dimensional network Q. For B = 1, 



the particles only interact with the external field; in this 
case only those particles that initially share at least one 
component of their vector states with the components of 
AP will converge to the field state MA The fraction of 
particles that do not share any component with M' is 
given by (1 — l/q) F ; thus the fraction of those particles 
that converge to M* is 1 — (1 — l/q) F . Figure 1 a shows 
both the numerically calculated values of S as well as the 
analytical curve of Sm versus q, for fixed B = 1. Both 
quantities agree very well, indicating that the largest do- 
main in the system possesses a vector state equal to that 
of the external field when B = 1. 

For intermediate values of B, the spontaneous order 
emerging in the system for parameter values q < q c due 
to the particle-particle interactions competes with the or- 
der being imposed by the field. This competition is man- 
ifested in the behavior of the order parameter S which 
displays a sharp local minimum at a value q*(B) < q c 
that depends on B, while the value of q c is found to be 
independent of the intensity B, as shown in Fig. 1 a. To 
understand the nature of this minimum, we plot in Fig. 1 
b the quantity a = S — Sm, as a function of q, where Sm 
is the normalized average size of the largest domain dis- 
playing the state of the field M f . For q < q*(B) the 
largest domain corresponds to the state of the external 
field, S = Sm, and thus a = 0. For q > q*(B), the 
largest domain no longer corresponds to the state of the 
external field but to other state non-interacting with 
the external field, i.e., S > Sm, and a > 0. The value 
of q*(B) can be estimated for the limiting case B — » 1, 
for which Sm ~ 1 — (1 — ^/l) F and the largest domain 
different from the field is S w 1 — Sm- Therefore the 
condition S = S M yields q*(B 1) = [l - (1/2) 1 /- F ] ~\ 
For F = 10 it gives q*(B — > 1) = 15 in good agree- 
ment with the numerical results. The order parameter a 
reaches a maximum at some value of q between q* and 
q c above which order decreases in the system and both 
5^0, Sm — ► 0. As a consequence, a starts to decrease. 

The collective behavior of the vector model on a fully 
connected network subject to an external field can be 
characterized by three phases on the space of parameters 
(q,B), as shown in Fig. 2: (I) an ordered phase induced 
by the field for q < q* , for which a — and S = Sm ^ 0; 
(II) an ordered phase in a state different from that of 
the field for q* < q < q c , for which a increases and 
S > Sm', and (III) a disordered phase for q > q c , for 
which a decreases and S — ► 0, Sm —> 0. 

For parameter values q < q c for which the system or- 
ders due to the interactions among the particles, a suffi- 
ciently weak external field is able to impose its state to 
the entire system (phase I). However, if the probability 
of interaction with the field B exceeds a threshold value, 
the system spontaneously orders in a state different from 
that of the field (phase II) . 

Continuous states based on bounded interactions pro- 
vide other instances of a noncquilibrium systems where 
induced and spontaneous order compete in the presence 
of an external field. Consider, for example, the bounded 
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FIG. 2: Phase space on the plane (q, B) for the vector model 
on a fully connected network subject to an external field, with 
fixed F — 10. Regions where the phases I, II, and III occur 
are indicated. 

confidence model 0] • It consists of a population of N par- 
ticles where the state of particle i is given by a real num- 
ber Ci <E [0, 1]. We introduce an external field M £ [0, 1] 
that can interact with any of the particles in the sys- 
tem. The strength of the field is again described by a 
parameter B £ [0, 1] that measures the probability for 
the particle-field interactions, as in the vector model. 

We start from a uniform, random initial distribution of 
the states of the particles. At each time step, an particle 
i is randomly chosen; 

1. with probability B, particle i interacts with the 
field M: if \d - M\ < d, then 

Cl +l = \{M + Cl); (1) 

2. otherwise, a nearest neighbor j is selected at ran- 
dom: if \Ci — Cj \ < d then: 

Ct +1 =C* +1 = ±(Ct + Cf). (2) 

The parameter d defines a threshold distance for interac- 
tion and the remainder we fix M = 1. 

We calculate the normalized average size of the largest 
domain S in the system as a function of 1 — d, for different 
values of B, as shown in Fig. 3(a). For B = 0, the system 
spontaneously reaches a homogeneous state Ci = 0.5, V«, 
characterized by S = 1, for values 1 — d < 1 — d c w 0.77, 
with d c w 0.23 [3]; while for 1— d > l — d c several domains 
are formed yielding S < 1. 

For B = 1 particles only interact with the field; in this 
case the value of M is imposed on the largest domain 
whose normalized size increases with the threshold, i.e. 
S = d. For intermediate values of B 7 the spontaneous 
order emerging in the system for values of 1 — d < 1 — d c 




FIG. 3: (a) S versus 1 — d for the continuous model for 
B — (stars); B — 0.5 (diamonds); B — 0.8 (squares); B = 1 
(circles), (b) a vs. d for B = 0.8 (circles) and B = 0.1 
(squares). The values of 1 — d c and 1 — d* are indicated for 
B = 0.8. Size of the system is N = 2500. 



due to the interactions between the particles competes 
with the order being induced by the field. The quantity 
S exhibits a sharp local minimum at a value 1 — d = 
1 — d* < 1 — d c , as shown in Fig. 3(a). In Fig. 3(b) 
we plot the order parameter a = S — Sm as a function 
of 1 — d, for different values of B. For 1 — d < 1 — d* 
the largest domain reaches a state equal to M, that is 
S = Sm, and thus a = 0. At 1 — d = 1 — d* , the state of 
the field no longer corresponds to the largest domain, i.e., 
S > Sm, and a starts to increase as 1 — d increases. For a 
small value of B, the quantity a reaches a maximum close 
to one, indicating that the spontaneously formed largest 
domain almost occupies the entire system, i.e., the field 
is too weak to compete with the attracting homogeneous 
state Ci = 0.5, Vi. However, when B is increased, the 
maximum of a is about 0.5, i.e., the attraction of the 
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FIG. 4: Phase space on the plane (1 — d, B) for the scalar 
model on a fully connected network subject to an external 
field. Regions where the phases I, II, and III occur are indi- 
cated. The dashed line in phase II separates region where the 
maximum of a — > 1 (below this line) from the region where 
a < 0.5 (above this line). 



field M = 1 increases and the size of the domain with 
a state equal to M is not negligible in relation to the 
size of the largest domain. In contrast, in the vector 
model the maximum a — ► 1 in the region q* < q < q Cl 
independently of the value of B. 

The value of d* in the scalar model depends on B and 
it can be estimated for B —> 1. In this case, Sm ~ d and 
S ~ 1 — d; thus the condition S = Sm yields d* « 0.5 
when B —> I. The quantity a reaches a maximum at the 
value 1 — d w 1 — c? c , above which disorder increases in the 
system, and both 5* and Sm decrease. As a consequence, 
a decreases for 1 — d > 1 — d c . 

As in the vector model, the collective behavior exhib- 
ited by the scalar model on a fully connected network 
subject to an external field can be characterized by three 
phases: (I) an ordered phase parallel to the field for 
1-d < 1 - d*, for which a = and S = Sm i= 0; 
(II) a ordered phase for 1 — q* < 1 — d < 1 — d c , for 
which a increases and S > Sm', and (III) a disordered 
phase for 1 — d > 1 — d c , for which tr decreases and both 
S and Sm decrease. Figure 4 shows the phase diagram 
on the plane (1 — d, B) for the scalar model subject to an 
external field. The continuous curve separating phases I 
and II gives the dependence d*(B). 

Short range interactions. To analyze the role of the 
connectivity on the emergence of an ordered phase or- 
thogonal to the external field, we consider a small-world 
network [To| , where the rewiring probability can be var- 
ied in order to introduce long-range interactions between 
the particles. We start from a two-dimensional lattice 
sites with nearest-neighbor interactions. Each connec- 
tion is rewired at random with probability p. The value 
p = corresponds to a regular network, while p = 1 
corresponds to a random network with (k) = 4. 




FIG. 5: S versus q in the vector model on a small world 
network with (k) = 4, B — 0.5, F = 3, for different values of 
the probability p: p — (empty circles), p = 0.005 (squares), 
p = 0.05 (diamonds), p = 0.1 (triangles), p — 1 (solid circles). 
Inset: S vs. p for fixed values q — 40 > q* and B — 0.5. 



Figure 5 shows the order parameter S as a function of 
q in the vector model defined on this network for different 
values of the rewiring probability p and for a fixed value 
of the intensity of the field B. The critical value q c where 
the order-disorder transition takes place increases with p, 
which is compatible with the large value of q c observed 
in a fully connected network. When the long-range in- 
teractions between particles are not present, i.e. p = 0, 
the external field is able to impose its state to the entire 
system for q < q c . Spontaneous ordering different from 
the state of the external field appears as the probabil- 
ity of having long-range interactions increases. The size 
of this alternative largest domain increases with p, but 
it does not grow enough to cover the entire system (see 
inset in Fig. 5). Increasing the rewiring probability in 
the scalar model also produces a behavior similar to the 
vector model. Thus, in systems whose dynamics is based 
on a bound for interaction, the presence of long-range 
connections facilitates the emergence of spontaneous or- 
dering not associated to the state of an applied external 
field. 

In summary, we have studied the collective behavior of 
nonequilibrium systems with non-interacting states and 
subject to an external field. We have considered two 
models on a fully connected network that share a com- 
mon feature: the existence of non-interacting states. In 
both cases we have found three phases depending on pa- 
rameter values: two ordered phases, one having a state 
equal to the external field, an another ordered phase, 
consisting of a large domain with a state orthogonal to 
the field; and a disordered phase. The occurrence of an 
ordered phase with a state orthogonal to the field is en- 
hanced by the presence of long range connections in the 
underlying network. We have verified that this alterna- 
tive ordered phase also appears when the models consid- 
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ercd here arc defined on a scale-free network. 

The emergence of an ordered phase with a state dif- 
ferent from that of an external field may be relevant in 
social systems as well as in many biological systems hav- 
ing motile elements, such as swarms, fish schools, and 
bird flocks [3| , whose dynamics usually possess a bound 
condition for interaction. Thus one may expect that this 
phenomenon should arise in large class of noncquilibrium 



systems in the presence of an external source for interac- 
tion. 
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